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(N 1 Abstract 



We determine the quark damping rates in the context of next-todeading order hard-thermaldoop 
summed perturbation of high-temperature QCD where weak coupling is assumed. The quarks are 



ultrasoft. Three types of divergent behavior are encountered: infrared, light-cone and at specific 



points determined by the gluon energies. The infrared divergence persists and is logarithmic 
whereas the two others are circumvented. 
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I. INTRODUCTION 



RHIC results seem to indicate that in the temperature range between one to two 
T c ~ 200MeV, the critical temperature of deconfmement from conventional hadronic states, 
there is an intermediary phase in which the deconfined quarks and gluons remain strongly 
interacting 111], a scenario introduced in WWW Since in this temperature regime the quark 
binding energy is about 4 GeV, one order of magnitude higher than the temperature itself, 
perturbative treatment based on quasi-free quarks and gluons is not adequate. Rather, a 
hydrodynamic description of a near-perfect liquid seems to work and the picture emerging 
is that of non-conventional bound states of quarks and gluons: diquarks qq, baryons qqq, 
three-gluon bound states ggg and polymeric chains qg . . . gq. A model of a strongly inter- 
acting classical chromo-electric plasma can account of much of the phenomenological results 

0. 

But the same quark binding energy decreases as the temperature increases away from T c 
Q], and here, starting from about 3T C , weak coupling between individual constituents may 
be expected. This is also supported by the behavior of particle susceptibilities: baryons 
and other bound states contribute in the intermediary temperature range, but only quarks 
survive in the high-temperature limit jsj]. 

Therefore, as the temperature increases, the picture emerging is that of a hadronic phase 
with very short inter-quark binding lengths changing at T c into the so-called strongly coupled 
quark-gluon plasma (sQGP), a near-perfect liquid phase of quark and gluon bound states 
with longer inter-quasiparticle binding lengths. As the temperature increases further, the 
sQGP changes into a weakly coupled quark-gluon plasma (wQGP) of quasi quarks and 
gluons with chromo- neutralizing isotropic Debye clouds Q]. 

If perturbative QCD is to apply in the wQGP phase, the problem becomes: how to 
organize it? It is known for some time that the standard loop-expansion would break at 
some order, depending on the quantity of interest under consideration j^, Q, [|. It has 
also proved inadequate when describing slow-moving particles since it does not reflect an 
expansion in powers of the coupling Regarding this last difficulty, an 

important improvement has been the dressing of the lowest-order propagators and vertices 
with the so-called hard thermal loops (HTL) 

0, 0, 0, 



19 



20, 



2l|. Since then, 



work has flourished and many phenomenological aspects of the presumed wQGP as well as 
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more theoretical aspects of high-temperature QCD have been addressed in the context of 

n 

lowest and next-to- leading orders of the so-called HTL-perturbation theory 22]. But HTL- 
dressed perturbative QCD is not itself safe from chronic problems, one important one being 
the non-screening of static chromo-magnetic fields at lowest order which plagues the theory 
with infrared divergences in next-to-leading order calculations. The determination of the 



chromo-magnetic correlation length may not even be accessible perturbatively 



221. 



The context is next-to-leading 



The present work finishes a calculation started in 
order HTL-dressed perturbative QCD and weak coupling is assumed. The aim is to complete 
the determination of the quark damping rates to second order in the ultrasoft external 
momentum. The coefficients of zeroth order are the ones already found in (24, 25]; they are 
finite and positive. The coefficients of first order we find are also finite and free from any 
divergence, but the coefficients of second order are logarithmically divergent in the infrared. 
This work is meant to be an additional contribution investigating the analytic infrared 
behavior of QCD at high-temperature. It is certainly more pressing nowadays to have a 
better theoretical understanding of the sQGP phase of hadronic matter, but ultimately, two 
issues have to be addressed: (i) Must we look for a quantum field theory description of the 
different phases of quarks and gluons? (ii) Do we want to reproduce the eventual phase 
transitions in this context? 

This article is organized as follows. After this introduction, section two recalls the es- 
sential results of |23j to which the present work is a follow-up. The analytic expression 
of the damping rates is given in the form of integrals over functions involving products of 
spectral distributions and their first and second derivatives. Section three takes up from 
these expressions and perform the integrals. The steps of the calculations are detailed using 
a generic form. The occurrence of divergences is discussed. Three types occur: infrared, 
light-cone, and divergences at specific points determined by the gluon energies. The infrared 
divergences are extracted and will stay, whereas the two other kinds are dealt away with. 



Section four summarizes the work and finishes the article with concluding remarks. 

II. QUARK DAMPING RATES IN HTL- SUMMED PERTURBATION 

We consider a theory with N c colors and Nf flavors. Imaginary-time formalism is used 
throughout. The euclidean momentum of the quark is P M = (po, p) such that P 2 = p\ + p 2 
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with po = (2n + l)vrT, a fermionic Matsubara frequency; n is an integer. Once all interme- 
diary steps are carried out, real-time amplitudes are obtained via the analytic continuation 
p = —iuj + + where uj is the energy of the quark. The temperature T determines the hard 
scale, gT the soft scale where g is the (weak) coupling constant, and g 2 T the ultrasoft scale. 
An infrared cutoff rj ~ g 2 T is introduced. 



A. Dressing the propagators and vertices 



The HTL-dressed quark propagator can be written as: 

*A F (P) = - [ 1+p A + (P) + j_ p A (P)] . (2.1) 

7^ are the euclidean Dirac matrices, 7 ±p = (7 ± 27. p) /2 and A ± = (P =|= P s ) _1 with: 

Wo ' 



Po (P) = ipo ~ ^Qo ( — ) ; D s (P) = p + ^ 
p \p J p 



m 



P VP 



(2.2) 



m f = V^f/^gT is the lowest-order quark thermal mass with Cf = (N 2 — 1)/2N C , and 
Qo {x 



X + 1 

I In -. The poles of A±(— iu, p) determine the dispersion laws u)±(p) to lowest 

order in g. The + sign is for real quarks and the — sign for the so-called 'plasminos' j^j], 
thermally excited quasiparticles. For soft quarks and using the notation p = p/nif, one has: 



u±(p) = m f 



1 1 o 16 , 1 . 32 

1 ± -p + ~p =F p H p ± 

y 3 F ^ 135 1 54 1 2835 



139 
12150 



p 6 ± O (p 7 ) 



(2.3) 



The quark damping rates 7 ± (p) are obtained by including in the dispersion relations 
the HTL-dressed one-loop-order quark self-energy *S(P). The inverse quark propagator 
becomes: 

A F X (P) = *A F X (P) - *£(P) . (2.4) 
The decomposition *£ = 7 *P + ^7-P implies: 



A F X (P) = - [ 7 ° (Po + * A)) + *7-P [D a + *D, 



(2.5) 



With the definitions 7±(p) = — Imf2± (p) with Q± the poles of Ap(— iO, p), and since *£ is 
(yf-times smaller than *A^\ we have to order g 2 T: 



7± (p) 



(2.6) 



w=u;± (p)+iO+ 
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where f± = D ^fD s , *f± = *Po = F *D S and stands for d/duj. Expanding the denominator 

in the above relation in powers of p using the expressions in (|2.2|l . one obtains: 

1 
2 



7± (p) 



l±^-^ 2 ±0(p 3 ) 



(2.7) 



aj=w±(p)+iO+ 



Therefore, determining 7 ± (p) amounts to calculating Im *£ (P) 



In HTL-summed perturbation 



14. 
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2l|, the one-loop-order quark self-enegy writes: 
S (P) = *£i (P) + *S 2 (P) . (2.8) 



The first contribution is: 



% (P) = -s 2 CVTr soft [T" (P, -Q; -If) *A F (Q) T" (-P, Q; if) *A j[i „ (if)] 



(2.9) 



a QED-like loop formed by two quark-gluon vertices connected by one gluon propagator and 
one quark propagator. The second contribution is: 



% (P) 



g 2 C f Tr soit (P, -P; if, -if) *A /Lti/ (if) 



(2.10) 



a purely hard-thermal-loop two-gluon-quark-antiquark vertex with one gluon propagator. In 
the above two expressions, if is the soft loop momentum, Q = P — if and Tr = J 



h) 



d 3 k 



with ko = 2nicT when bosonic or kg = (2n + 1) nT when fermionic. The subscript 'soft' 
means that only soft values of if are allowed in the integrals; hard values have already 
dressed the propagators and vertices. 

To be complete, we give the expressions of the gluon propagator and the vertices involved 
in flUD and (ETTUl) . In the strict Coulomb gauge, *A 00 (if) = *A, (if), *A 0l (if) = and 
*A tj (if) = (Sij - kkjj *A t (if) with *A, and *A t given by: 

1 . 1 



"A^ (if) 



k 2 - Silt (if) 



*A, (if) 



(2.11) 



where SU^K) = 3m 2 Qi(^f) and 5U t (K) 



if 2 - 5U t (if) ' 
[Q 3 (f ) - Qi(^) - |] are the gluonic 



hard thermal loops. Here Qi(jr) is a Legendre function of the second kind and m g = 
yiVc + Nf/2gT/3 is the gluon thermal mass. Finally, the HTL-dressed vertices *T are as 
follows: 



T^(P,Q;P) = Y + m} 
*f " U (P, -P; if , -if ) = -2m 2 , 



4tt PSQS 1 



f J 4tt PS (P + K)S (P - if)5 
In both (j2.12|) and ()2.13|) . S = («, s) and fi s is the solid angle of s. 



(2.12) 
(2.13) 



B. Analytic expressions 

Let us recall the main results of [23]. Henceforth, the quark thermal mass m/ is set 
to one. Since there remains another soft mass in the problem, m g , we define the ratio 
m (JV C , Nf) = m g /mf = |y I N ^ N ^+^/ 2 ) ^ j s eaS y ^ Q see we always have m > 1. Since 
the gluon propagator is taken in the strict Coulomb gauge, there are uncoupled longitudinal 
and transverse contributions to *£i and *£ 2 . There is a further split in *Ej due to positive 
and negative helicity quarks. Six contributions in all, each being treated separately. 

After a preliminary manipulation of the gamma matrices, we perform the expansion 
'ps = ip^ 1 — |% ~ ^ (^° 3 ) * n or der to integrate analytically over the solid angle Q s . 

This expansion is valid in the region p < \ipo\, a condition always satisfied before analytic 
continuation because po = (2n + 1) 7rT and p ~ g 2 T, and after since for ultrasoft momenta, 
ip = rrif + 0(p) ~ gT, see f!2.3|) . A subsequent expansion in powers of p of functions of 
q = |p — k| is necessary in order to carry out analytically the integration over the internal 
three-momentum solid angle Qk- In the process of these angular integrations, special care 
must be taken when rotating the gamma matrices 



Next we perform the Matsubara sums. For this, the spectral 



decompositions of the 



dressed propagators and Q (K) with k fermionic are used J2l|,|27|, 28]: 

A £ {k ,k) = / dre ikoT dup £ {uj,k) (1 - fi{uj)) e- UT ■ 

JO J~oo 
fl/T p+oo 

A t (k ,k) = / dre ik0T dup i {uj,k)(l + n{uj))e~ u}T ; 

J0 J -oo 

Q (ik /k) = / dre ikQT / du p {uj, k) (1 - n(uj)) e~ WT . (2.14) 



JO J-oo 

e stands for + or — and i for / or t. The functions n(uj) and n{uj) are the Bose-Einstein 
and Fermi-Dirac distributions respectively, and the rho's are the spectral densities, to be 
displayed in the next section. Before implementing the spectral decomposition, it is first nec- 
essary to rearrange terms in such a way that products of at most two functions necessitating 
a spectral decomposition occur to ensure the obtainment of only one energy denominator 
just before the extraction of the imaginary part. The steps of the calculation must also 
ensure the Matsubara frequency iko appears only in the numerator of fractions. Hence, 
using (j2.14|) . the sum over k Q can be performed, yielding a delta function that automatically 
removes one integration over one imaginary time while the other integration produces an 



energy denominator. At this stage, every ipo is to be replaced with (2n + 1) ttT except in 
the energy denominator. 

Now the analytic continuation to real energies ipo — > uJ±(p) + i0 + can be taken. The 
extraction of the imaginary part becomes straightforward using the relation 1/ (x + i0 + ) = 
Pr (l/x) — i7rS(x) where Pr stands for the principal part. Further rearrangements are made 
and, according to the definitions in (|2.7jl . the quark damping rates are given by the following 
expressions: 



7± (p) 



9 2 CfT 



in 



a G ±a 1 - + a 2 — + 0(p i ) 



(2.15) 



Recall that we have set to/ = 1. The coefficients aj are given by the expressions: 



Of) 



dk 



+oo 



du 



oo 
+oo 



duo 



UJ 



7 /o5: 



poo p+oo p+oo jl^I 

ai= dk duj —[fi-fo dj\ 5, 

Jrj J —oo J —oo 

poo r+oo r+oo j i 

a 2 = dk du — [f 2 - fid u - f (3 d u - dl)] 5, (2.16) 

J r/ J —oo J —oo " t/ 

with 5 = 5 (1 — uj — uj'). The three functions /j = /j (a;, a/; fc) are given by the following 
rather long expressions j^j: 



1 



(i _ £fc + W )2 pep ; + - (l + 2e:/c + /c 2 - ^ 2 ) p eP ; 



+ - (k 2 - uj") p p' t . (2.17) 



This expression is the one obtained in 



24, 



251 ] for the non-moving quark damping rates. In 



the above expression and the two subsequent ones, the notation is as follows: p £ stands 
for p £ q{uj, k) and p' lt for p lt (uj', k). Remember also that d w stands for the partial derivative 
d/duj. The next function f\ (uj,uj'; k) is given by: 

2s 



2A; 2 (-1 + k 2 - 2eku + uj 2 ) p e p\ + 



3 + 2ek + 4k z - k 



PePt 



- (2 + Ask + 2A; 2 ) uj + + 4 + 2ek + 2fc 2 ^ tu 2 + 2cj 3 - 0£ + lj , 
+efc 2 (l- efc +w) 2 p £ <9 fe p; + [| + 2A; + 3ek 2 + 2k 3 + £ -k" -{e + 2k + ek 2 ) uj 2 + £ -uj A ^p £ d k p' t 
-1 [k 2 -uj 2 + 2^)p p' t -2k 2 e (uj) 6 {uj 2 - k 2 )p\ +^(uj 2 - k 2 )p d k p' t + 2ujp d k p' l , (2.18) 

where d k is a short notation for the partial derivative d/dk and e (uj) is the sign function. 
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The last function / 2 (no, uj'; k) to display is the longest of all three and is given by: 

f2 = ~ k2 - 6ek3 - \ kA - ( 6ek - 6k2 + 2ek3 ) u + ( 9 + fc2 ) ^ 2 + 6ekuj3 - 



e=± 
XPeP? + 



+ 



9 

6 2e 
k? + ~k 



Ue 8 



19 



fc 2 -6efc 3 + fc 4 + 



-3 25e 



fc 2 



+ 



•10+6efc + 9fc 2 -3efc 3 w 



23 - 6ek + k 2 )uj 2 + 



6 22e 



A: 2 



k 



-6 + 6ek)u 6 + 



3 

2P 



12e 



-5 a; 



(3 3 "\ 3 "\ 

- + -£)a; 8 + t^ 6 J P e P^ - fc (9 - 14eA; + 5A; 2 + (12 - 2ek - 6k 2 ) uj - (3 - 12eJfe) a> 2 



> 2 + fc , 

6a> 3 ) p e 9 fc p; + ( ^ + 4 £ + 7k + 8ek 2 + h 3 - 



+ 6e + 6k + 6efc 2 + 3fc 3 a; 



-(| + As + 5k\ uj 2 +{^ + 6s + + - |a; 5 ) p^pj - ^ 2 (1 - e* + P £ <9 2 P ; 



> 2 



)PoPi 



+ g + 3e* + \k 2 + 3ek 3 + - I (1 + 2eJfc + A; 2 ) c 2 + p £ d 2 k Pt 

+ ("3 + Qku + j^uj 2 - -a; 3 ) p d k p' t + Skp^tpl - A [k 2 - c 2 ) p 9 2 P ; 

+ 12A; 2 e (w) 5 (cj 2 - A; 2 ) p\ - 6uje (w) 6 (uj 2 - k 2 ) p' t - 6k 2 ue (u) 3^5 (uj 2 - k 2 ) p\ . (2.19) 

Note that, since only soft values of u and uj' are to contribute, we have made use of the two 
approximations n(uj) ~ 1/2 and n(uj) ~ T/uj. 

To complete the calculation of the damping rates, it remains to perform the integrals 
over the frequencies uj and uj' and then over the momentum k. These integrations are not 
straightforward and necessitate numerical work. But before that, one has to hunt down 
divergences, particularly the infrared ones. Also, the dimensionless parameter m(N c , Nf) is 
implicitly present in the spectral densities p l t and so, each case (N c , Nf) has to be treated 
separately. All this is discussed in the next section. 



III. INTEGRATION AND DIVERGENCES 

In this work, we will consider three cases regarding the number of colors and flavors, 
namely, (N c , Nf) = (2, 1), (3, 2) and (3, 3). All the terms in the expressions of the coefficients 
do, ai and a<i in (|2.15jl have the following generic structure: 

/•OO /■+OO /-00 J / 

I rs = dk duj —g(uj,k)p(uj,k)d r k p f (uj' } k)d s J(l-uj-uj'), (3.1) 

J T] J— OO J— OO 



where g (uj, k) is a polynomial in uj with coefficients functions of k. The density p (uj, k) 
stands for the quark spectral functions given in ()3.6|) below or for the spectral distributions 
given in (|3.9|) . The density p' (uj', k) stands for the gluonic spectral functions p l t (uj' , k) given 
in (J3.2j) . The two indices r and s designate the rth and sth derivatives with respect to k 
and uj respectively and are such that r, s = 0, 1 or 2 with the condition r + s < 2. 

First we give the expressions of the spectral functions. The gluonic spectral densities are 
as follows: 



p lt (uj, k) = i ut {k) [8 (uj - uJi tt (k)) -6(u + vi,t(k))] + (3 lt (u, k)6(k- . 



(3.2) 



The gluon energies ooi t t(k) are the poles of the longitudinal and transverse gluon propagators 
*Aij (K) given in (|2.11|) . The longitudinal and transverse residue functions are given by: 
utWitf-ulik)) n . uj t (k) (ujf (k) - k 2 ) 



' k 2 (3m 2 - lu 2 (k) + k 2 ) ' 
The longitudinal cut function is given by: 



3m 2 LU 2 (k) - (tu 2 t (jfc) - k 2 f 



-3m 2 w 



2k 



and the transverse one by: 
P t (k,u) = - 



(3m 2 + F-3m 2 fln(^)) 2 +(^)- 
3m 2 uj(k 2 — uj 2 ) 



(3.3) 



(3.4) 



2k 3 



(^-k 2 -^ 2 (^ + ^p^ln(|^))) 2 + 



3Trm 2 uj(k 2 —U! 2 ) 

ifc 3 



(3.5) 



The quark spectral functions are: 



p± (uj, k) = t±(k)6 (u - u±(k)) + l T {k)5 (uj + uj ± (k)) + (3 ± (uj, k)0(k- \u\) 
where the quark residue functions are: 

l±(k) = -\ K (k) - k 2 ) , 



and the quark cut functions given by: 
P± (w,A;) = - 



—(k =F uj) 



2k 2 



2k 2 



(3.6) 



(3.7) 



(3.8) 



(uJTk + ^((k T uj)\n(^) + 2k)) 2 +[ 

The other spectral distributions appearing in the functions /o, fi and f'2 are as follows: 

p (uj, k) = ^9 (k 2 — uj 2 ) ; p 1 (uj, k) = e (uj) 5 (k 2 — uj 2 ) ; p 2 (uj, k) = e (uj) 8^28 (k 2 — uj 2 ) . 

(3.9) 



To show how to carry with the generic integral I rs , it is appropriate to chose the spectral 
function p (co, k) as a quark spectral density. This example is typical and general enough 
to encompasses all the major difficulties and subtleties we encounter throughout the work. 
Using the expressions in (|3.6|) and (J3.2j) . we see that there are four kinds of contributions: 
pole-pole (pp), pole-cut (pc), cut-pole (cp) and cut-cut (cc). Let us start with the pole-pole 
contribution, which typically writes as: 

poo poo poo J I 

III = dk du —g (u, k) i ±£ (k) 5 (u =F co £ ) d r h [ 3m (k) 5 (co' ± u h t)W (1 - o> - «/) ■ 

J rj J — oo J — oo 

(3.10) 

The integrand is nonzero only at the intersections of the supports of the delta functions, 
namely, when u ± co £ = 0, co' ± co^ t = and u + to' — 1. This is satisfied at the points k 
for which we have ±co £ (k) = 1 ^co^t (k). But for the three cases (N c , Nf) considered in this 
work, these curves never intersect. Hence we can write: 

III = 0- (3.11) 
Next we look at the term pole-cut. It is of the form: 

poo poo poo J J 

III = dk du —g (co, k) h S (co ± co £ ) dl r) [p lit (co', k) 6 (k - \co'\)] d^5 (1 - co - u') . 



(3.12) 

With the aim of making the illustration clear, let us consider first the specific case 
r = 1 and s — 0. The integration over co' becomes trivial. The derivation with re- 
spect to k yields two terms: g (co, k) i £ (k) d^(5 lt (1 — u, k) 5 (co ± co £ ) 9 (k — |1 — co\) and 
g (co, k) 3 e (k) /3 lt (1 — co, k) 5 (co ± co £ ) 5 (k — \ \ — co\). The second term is always zero be- 
cause it imposes the conditions co £ (k) = ± (1 =F k), which are never satisfied kinematically. 
In the first term, the integration over co is straightforward and the theta function imposes 
the constraints —k < 1 ± co £ (k) < k on the momentum integration. Only the minus sign 
can be satisfied, and for all values of k. Hence we write for this contribution: 

j " = r dk r^wf ^ ^ ^ k) .« • (3 - i3) 

The case r = 2 and s = is carried out in a similar way. Because of the second derivative in k, 
there will be three contributions: (i) g (co, k) % £ (k) d\{5 { t (1 — co, k) 5 (co ± co £ ) 6 (k — |1 — co\). 
The integration over co is trivial and the theta function imposes the minus sign with 
no constraints on k. (ii) 2g (co, k) } £ (k) dk(3 l t (1 — co, k) 5 (co ± co £ ) 5 (k — |1 — co\) and (iii) 

10 



g (uj, k) i £ (k) (5 l t (1 — uj, k) 5 (uj ± uj e ) dkS (k — \1 — ou\). These two contributions cancel be- 
cause of the same kinematics as in the previous case. Then we can write: 



! ™ = l <tt r^ fc( * ) ^ (w -* ) U~e.r (3 ' 14) 

The other cases for r and s are worked out in similar steps. When a derivative with respect 
to uj intervenes over 5 (1 — uj — uj'), it is first transformed into an integral over uj' and then 
brought back onto the other functions. The rest of the steps are straightforward. We can 
write generically: 



poo 

i£ =(-)'/ dkg(uj £ ,k) d£ (k)dZ 



(3.15) 

LO = l—L0 e (k) 



UJ 

The notation d r k ,s means derive r-times with respect to k and s-times with respect to uj 



A. Infrared behavior 



The generic integral 1*% is present in all three coefficients ao, ai and a 2 . From the explicit 
expression of i e (k) given in (j3.7J) and those of P tt (uj, k) given in ()3.4|) and (j3.5|) . one carries 
out the fc-integration in f)3.15j) for the different specific functions g (uj, k) intervening in the 
different contributions to fo, fi and fi given in relations (j2.17j) . ()2.18|) and ()2.19|) respectively. 
Most of the work is numerical. For the coefficients ao and a±, no particular difficulty arises 
and the limit rj — > is safe. However, the coefficient a 2 requires special attention since the 
infrared limit is sensitive. To see this concretely, take the explicit example of p + d1p[ in 
f 2 with g (uj, k) = —3k 2 (l — k + uj) 2 . Let us write the integral as I 2 q = dkF (k) with 
the integrand F (k) = -3k 2 (l — k + oJ+fl+l (1 - uj + ) dlf3 l (uj, k)\ u=1 . With the explicit 
expressions of uj + (k) given in (|2.3jl . j + (k) in ()3.7|) and (3 l (uj, k) in ()3.4|) . we can perform a 
small-fc expansion of the integrand F (k) to find for the different cases of N c and Nf. 

4486 

(N c ; N f ) = (2, 1) — > F (k) = — — — + 0.4486 + 0.2248A; + O (k 2 ) ; 

K 

(N c ; N f ) = (3, 2) — > F (k) = — — — + 0.5428 + 0.2466A; + O (k 2 ) ; 

rC 

4431 

(N c ; N f ) = (3, 3) — ► F (k) = — — — + 0.4825 + 0.2223A; + O (k 2 ) . (3.16) 

K 

The 1/k behavior indicates a logarithmic infrared divergence. This latter is extracted in the 
following manner. Writing F (k) = a/k + finite, we have: 

J 2 P C = -alnr/ + alnA + ^ dk (f (k) - + dkF (k) . (3.17) 
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In short, we have split the original integral into two: one from i] to an arbitrary number 
A and another one from A to oo. The second integral is finite. From the integrand of the 
first one we have subtracted a/k, which makes the integral safe in the infrared and the limit 
i] — > can be taken. The term a/k has to be integrated by itself from rj to A. As said, the 
choice of A is arbitrary and we must (and do) check that the result in independent of it. For 
this specific example, we obtain the following results: 

(N c ; N f ) = (2, 1) — ► 1% = 0.4486 In 77 - 0.7641; 
(N c ; N f ) = (3, 2) — ► III = °- 4985 lnr 7 - 0.8557; 

(N c ; N f ) = (3, 3) — ► III = °- 4431 ln ?7 ~ 0.7539. (3.18) 
All other infrared-sensitive contributions are worked out in a similar manner. 



B. Other circumvented singularities 



The upcoming cut-pole contribution is not sensitive to r\ but has other difficulties we 
discuss now. Consider then: 



TCP 

rs 



-oo poo j 1 

<lk I duo / —g (u, k) /3 e (uj, k)0(k— \u\) d r k faS (00' - c*)] d s J (1 - u - J) , 



(3.19) 

where i stands for I or t. Let us be a little more specific and consider the case r = 2 and 
s = 0. The integration over uj' is trivial. We obtain: 

foe roo j 

1% = dk [g (cu, k) /3 e (w, k) 9 (k - M) y + 2^d k + M dl] 5(1 -u- uj,)} . 

Jr) J -00 

(3.20) 

Here, 3^ (k) is the first derivative and 3" (k) the second derivative of the residue function. 
The kinematics imposes 1 — k < u>i (k) < 1 + k, conditions satisfied for all k from the lower 
bound ki to 00, where k{ is the solution to the condition oj-i (k) — 1 + k. In the last two 
contributions, the /c-derivatives over the delta function are brought onto the other functions. 
With some algebra and using the fact that lim^-^ (3 £ (uj, k) = to eliminate some of the 
intermediary terms, we can write: 

POO 

I?o = / dk [fll - (itu'l + 23^) + uu'fdl] [g (oj, k) P e (oj, k) I (1 - u;)],,^ 



+ 



fau?g(l -Uj,k) 

OJi I 1 - u' ; I 



2d^ t (w.AOUi - 



P' e (l-Ui,k) 



(3.21) 
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In this expression, f3' e (1 — uji, k) indicates the total derivative with respect to k. Now we 
must be careful with this expression because the derivatives of the quark cut functions at 
ki are infinite. This is not related to the infrared limit. However, putting these infinities 
together cancels them. More specifically, if we call y = k — ki, then the singular behavior of 
the integrand around ki coming from the derivatives of the quark cut functions has the form 
A(\ny,y) /y 2 + B(\ny,y) jy, where A and B are complicated but computable functions. 
The task is then to systematically extract this singular behavior and put it together with 
a similar one coming from the terms not under the integral sign. It turns out that in each 
case and for every term, when put together, all the singularities cancel. 

C. Light- cone behavior 

Last we turn to the cut-cut contribution to I rs . Here a different kind of potential singu- 
larities arises, light-cone, but circumvented too. The generic term has the following form: 

pOO pOO pOO T / 

I c r c s = dk du —g(u,k)P e (u,k)8(k-\u\) 

J rj J —oo J —oo ^ 

xd r k k)6 (k - d™5 (1 - u - a/) . (3.22) 

Here too the case r = 2 and s = is typical enough. Performing the integral over uj' trivially, 
we have: 

POO poo J 

I2I = dk - g (uj, k) [3 £ (uj, k) 9 (k - 

JO J-oo i ~ UJ 

x [d 2 Ml -u,k)6(k-\l-u\) + 2dMl -u,k)S{k-\l-u\) 

+ Pi(l-u,k)d k 6 (k -\l-cu\)}. (3.23) 

The 00 term is constrained byl — k < uj < k with k > 0.5. We obtain the contribu- 
tion J °^ dk f*_ k dug (oo, k) (3 £ (uj, k) 9^/3^(1 — uj, k)/ (1 — u). The &5 contribution is left with 
one integral over k; it reads 2 J °^ ^-g (1 — k, k) f3 £ (1 — k, k) dkfi^l — uj, k)\ u=1 _ k . The Qd k 5 
contribution is treated similarly and yields — J^° 5 ^-dk [g (uj, k) f3 £ (uj, k) (3^1 — uj, k)] u=1 _ k . 
However, it is easy to see that we can get rid of the second and third contribu- 
tions by performing one integration by part in the first contribution. It will yield 
— J °^ dk f*_ k jz^d k [g (uj, k) (3 £ (uj, k)] d k P i (1 — uj, k) plus the opposite of the second and 
third contributions. Hence, changing the integration from over uj to over 1 — uj, we can 
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write: 



f°° f duo 

I20 = ~ dk —d k [g (1 - u, k) (3 £ (1 - u, k)] d k (3 t (u, k) . (3.24) 

</0.5 Jl-k ^ 

But extra care must be taken because lim^fc dkPi (u>, k) in infinite. This situation is 
handled as follows. Given the explicit expressions (|3.4|) and ()3.5|) of the longitudinal and 
transverse cut functions respectively, one can see that the divergence in (u, k) comes 
from l/( = 1/ (k — uS) and Y^ 1 = \n(. Then, writing J|q = J^dk J^ k 1 d(G (k, (;Y), we 
expand the integrand G (k, £; Y) in powers of £: 

G (k, C; Y) = G_! (*, F) /C + G (/c, Y) + d (A;, y) C + O {C) • (3.25) 

The singular term G-i (k,Y) /( is singled out in the following manner: 

r+oo r2k — l poo r0~ j-y 

I?o= dk dC[G(A ; ,C;l/lnC)-G_ 1 (fc,l/lnC)/C]+/ dk —G^(k,Y). 

JO.5 JO JO.5 Jln- 1 (2k-l) Y 

(3.26) 

In this expression, the first integral is finite. In the second integral, knowing that 
G-i(k,Y) = O (y 3 ) close to zero, there is no divergence anymore. Hence the whole i|g 
is finite. The same trick is used for all the cut-cut terms sensitive close to the light-cone and 
all singularities are circumvented. 

IV. RESULTS AND CONCLUSION 

All the terms in f , fi and fi given in ()2.17|) . (j2.18|) and ()2.19|) respectively have to be 
calculated. There are no additional subtleties to mention. Assembling all the partial results 
together and reintroducing the quark thermal mass m/ such that p = p/rrif and fj = r]/mf, 
we find the following final results: 

3^ 2 T 





N f ) 


= (2,1)- 


- > T± (p) 


(N c ; 


Nf) 


= (3,2)- 


7± (p) 


{N c ; 


Nf) 


= (3,3)- 


- > 7± (p) 



64tt 

12tt 
S 2 T 



[5.6978 =F 1.0452p— (6.0427 In 77 - 8.4684)p 2 +O (p 3 )] ; 
[5.6344 =p 0.9492p - (6.7141 In 77 - 7.7539)p 2 + O (p 3 )] ; 
. [5.7057 =Fl.0568p- (5.9681 In fj - 8.5536)p 2 + O (p 3 )l . 

127T 

(4.1) 

Remember that the present calculation of the quark damping rates is done in the context 
of next-to-leading order hard-thermal-loop summed perturbative QCD at high temperature 
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where weak coupling is assumed. It is intended to be an additional contribution to probe the 
analytic properties of finite-temperature QCD, particularly in the infrared. As mentioned 
in the introductory remarks, it has no direct relevance to the physics of the strongly coupled 
quark-gluon plasma, but it may be useful if, eventually, quantum chromodynamics is still 
believed to be the bedrock of all theoretical hadronic physics. 

One peculiarity of the results (j4.1j) is that the logarithmic sensitivity to the infrared is 
found only in the coefficients of p 2 ; the zeroth and first order coefficients are safe. The 
calculations have encountered other potential divergences: at the light-cone and at specific 
points determined by the gluon energies. But all these have been circumvented. 

The persistence of the infrared divergence is most probably attributed to the non- 
screening of the static chromo-magnetic fields at lowest order. In a separate work [29], we 
have considered the context of scalar quantum electrodynamics, where calculations are more 
fluid, and have explicitly shown that an early momentum expansion like the one we perform 
in [2^| and this work gives exactly the same results as a late expansion, after the Matsubara 
sum and analytic continuation to real energies are done. Remember that only soft values 
of the internal momentum have been included in the integration. It seems then that the 
ultrasoft region is important and needs to be investigated with probably non-perturbative 
means. 
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